Superconducting decay length in a ferromagnetic metal 
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The complex decay length ^ characterizing penetration of superconducting correlations into a 
ferromagnet due to the proximity effect is studied theoretically in the frame of the linearized Eilen- 
berger equations. The real part and imaginary part ^2 of the decay length are calculated as 
functions of exchange energy and the rates of ordinary, spin flip and spin orbit electronic scattering 
in a ferromagnet. The lengths ^1,2 determine the spatial scales of, respectively, decay and oscillation 
of a critical current in SFS Josephson junctions in the limit of large distance between supercon- 
ducting electrodes. The developed theory provides the criteria of applicability of the expressions for 
^1 and ^2 in the dirty and the clean limits which are commonly used in the analysis of SF hybrid 
structures. 
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The decay length ^ is an important material parame- 
ter which characterizes the scale of penetration of super- 
conducting correlation into a non-superconducting ma- 
terial across an interface with a superconductor. The 
critical current Ic in a Josephson junction scales expo- 
nentially with the distance between the superconducting 
electrodes L if i is larger that Ic oc exp{— L/^} . In 
nonmagnetic materials the decay length is a real num- 
berjwhile in ferromagnets ^ is a complex number (see 
0-01 for the reviews). In particular, if the condition of 
so-called dirty limit is fulfilled in the F metal, the decay 
length is 



{ttTY +H^±ttT 



(1) 



where Dp and H are the diffusive coefhcient and the 
exchange field in a ferromagnet, respectively. In the clean 
limit 

i^' = ^o'+i-\ ^0' - — , ^2' = Ch' = — , (2) 

Vp Up 

where vp is the Fermi velocity in a ferromagnet and t is 
the electron mean free path. From iQJ, (0) it is clearly 
seen that for dirty materials ^2 > ^1, and in the limit 
of large H » ttT, the characteristic lengths are nearly 
equal ^1 ~ ^2- In the clean limit these length scales ^1 
and ^2 are completely independent. 

The existing experimental data obtained up to now 
in SFS Josephson iunctionsp|-p^ can be separated into 
two groups depending on whether weak or strong ferro- 
magnet was used for junction fabrication. To be con- 
sidered as a weak ferromagnet, the dilute ferromagnetic 
alloys (e.g. Cui-^; Ni^,,) should be in the range of concen- 
tration close to the critical one {x « 0,5). The electron 
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mean free path in these alloys is very small providing the 
fulfillment of the dirty limit conditions. As a result, the 
observed relation between the decay (^1) and oscillation 
(^2) lengths ^2 ^ Cii is close to that following from Q- 
It is necessary to point out that in some experiments jl2l | 
the observed difference between ^2 and ^1 is so large that 
it can not be explained by temperature factor in |^ only 
and spin-dependent scattering processes should be taken 
into account (12.],,17J. 

Contrary to that, in the structures with strong 
ferromagnet [l6| (Ni, NiaAl), the relation between ^1 
and ^2 is just the opposite and large ratio ^1/^2 ~ 10 was 
observed in NisAl [1^1 . Therefore more complex model 
should be developed for the data interpretation. 

Most of previous theoretical work on SF hybrids was 
performed assuming the dirty limit (see 0-0), and only 
first order corrections to the decay length in small pa- 
rameter l^H ^ 1 were discussed in jl^ il9| . Properties 
of SF structures in the clean limit were also studied in a 
number of papers, see e.g. |^ |22, 113, 13 ■ The purpose 
of this work is to develop general theory describing the 
decay length ^ in a ferromagnet for any relation between 
^0, and i. 

To do this we consider a generic SFS Josephson junc- 
tion with arbitrary transparency of SF interfaces and 
large thickness of the F layer L ^ ^1. It is well known0- 
that the critical current of this structure should fall 
exponentially with L 

/c = /oexp{-i/^}. 

Here the prefactor Iq depends on physical properties of 
SF interfaces and the nearby S and F regions, while ^ 
depends onl y o n bulk parameters of F material and can 
be obtained|26|.[27l| as the solution of linearized quasi- 
classical Eilenberger equations f25j. These equations are 
valid at the distances from the interfaces larger then ^ 
and have the form HI], 0,0 



2 



a: 



^10 



{^o'±^G')f± + cose—u 



i^f'f ((/±> - f±) + Ceff ih - /±) . (3) 



soeff \JT 



(••■)sin 



(4) 



(5) 



Here 6 is the angle between the direction of electron ve- 
locity vf and the a;— axis, which is oriented perpendicu- 
lar to the interfaces, /± = f±{x, 0) are the quasiclassical 
Eilenberger functions describing the behavior of spin up 
and spin down electrons in the presence of exchange field 
H oriented parallel to the SF interfaces. The parame- 
ters = vpTso,£z = vpTz £x = vpT^, are the electron 
mean free paths for parallel and perpendicular to the di- 
rection of H magnetic scattering, while £so = vpTso is 
the electron mean free path for spin orbit interaction. 
Solution of Eq. © has the form 



/±(x, 0) - C±{9) exp , - + ^C: 



(6) 



Averaging in (|10|l over angle 6 we get the system 
of two equations for {C±{0)) . Its compatibility condition 
results in the equation for the effective decay length ^^f / 



tanh ■ 



4o ± 



(11) 



It is clearly seen that if the effective spin orbit interac- 
tion is so strong that £Joeff — ^h^' then the right hand 
side of Hll|l is real. Therefore in this case Eq. Ijlll) pro- 
vides us by two solutions for while ^2^^ = 0. It is 
necessary to mention that in the absence of ferromagnetic 
ordering {H — 0) due to degeneracy in spin orientation 
the critical current must not depend on £ soeff - In this 
situation only the root of equation corresponding to the 
'+' sign in Eq. should be considered 



tanh 



£~^ 

^eff 



Co 



^eff 



(12) 



which provides the largest value of the decay length. 
Solution of Eq. lfTT|) 



tanh 



"■eff 



Co 



£: 



ff 



2£ 



-1 

soeff 



(13) 



where ^ is the effective decay length independent on 9. 
Substitution of ^ into © provides the system of two 
equations for C±{6) 

(Co^^ + iCH')C+iO) - cos9C+i9) = 



= C// i{C+i9)) - C+i9)) + Ce// {C-{9) - C+{9)) (7) 

= C// {{C-{9)) - C-(0)) +Ce// {C+{9) - C^{9)) (8) 
Solution of these equations has the form 

{C+{9))K-J+£-,\jf {C^{9)) 



C+{9) = 



^eff 



(a-a- 



'^soeff 



(9) 



with the smaller ^ = ^12 also exists at finite H. (In the 
limit H ^ the prefactor before this exponential solu- 
tions goes to zero 17] providing independence of the crit- 
ical current on ^12). At £soeff = Ch these two lengths, 
are equal to each other, = ^12. With further H in- 
crease the right hand side of Ec ljllll becomes complex 
and Ea. ljlll) can be rewritten as 



tanh ■ 
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^^2 



-2 

soeff 



(14) 



The sign '-' in Ea. ljll|l simply provides the equation for 
the complex-conjugate solution of Ea. lfnji . 

In the limit £eff <C ^ one can expand the hyperbolic 
tangent in H14f) in series keeping three first terms and get 



sr. 




^<^ff 
^20^+ 



(15) 



C^{9) = 



{c^{9)) K'+Ceff {c+m 



^eff 



(a-a- 



^soeff 



ClO — Co + ^eff + ^ soeff 



(10) 



r+ = 



3r- L , 1^4// 
10 




- 1 



i + r, 



/2 
S20 



± 



^01 



(16) 
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The expressions in the square brackets in (|15|l . p6l) give 
first order corrections to the dirty Umit formula 17] for 
^1 and ^2- This approximation vaHd if 



soeff 



Iff 

(17) 

In the hmit ^o,^soe// ^ the expression ^ = 
(1 - |ii?r), T = l/vF, foUows from Eqs.lUni, (CHIl- 

This formula was obtained before in Ref. [l^ il9t i2Q] 
and can be interpreted as a complex correction to the 



diffusion coefficient, Dp 
In the clean limit 

A > max < In \/34, In 



Dp 1 



HHt) 




^soeff 



(18) 



A = l 



Co 



soeff 



in the first approximation we may put the hyperbolic 
tangent in l|14|l equal to unity and get 



-1 

10 ' 



^2^' = ^2^ (19) 

It is clearly seen that for iJoeff ~^ ^ ^^^^ formula trans- 
forms into Eq.lPll- In the next approximation it is easy 
to get that the corrections to 1)19(1 



cr' = cro'-2pcxp(- 
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eff 



I exp 



2e. 



eff 



;io 



(20) 



(21) 



where 



P = ^20^ sin 



q = sin 



21. 



eff 



,20 



2L 



eff 



C21 



20 



ClO^ COS 



^20^ COS 



21 



eff 



>20 



21 



eff 



20 



are oscillating functions of ^h- 

Eg. 1(14(1 is equivalent to the system of equations for ^1 
and ^2 



, 2er^ / 2^0 

coth — ^ a cos ^4- 



10 



"■eff 



^eff 



arctan ^ ) (22) 

croV 



From the structure of equations 1(22(1 . I(23() it follows 
that increase of £^20 leads to increase of ^2^^. This, in 
turn, results in increase of the second negative item in 
right hand side of 1(22(1 . Since ^j"^ must be a positive 
value increase of should be accompanied by a jump, 
at a certain point, to the positive branch of cos(a;) lead- 
ing to a discontinuity of S,2^{£,2()) dependence. This con- 
sideration is proved by numerical solution of p4(l (see 
Figs.1-3). 

Figures 1 and 2 show the dependencies of Cr^(C20^) 
and £,2^ {^20) ~ ^20 calculated for fixed values of param- 
eter ■ Open triangles and circles in the figures show 
the asymptotic dependencies (IT^ . l(TC|l and (|2DJ), l(?T|l . 
respectively. It is clearly seen that in the parameter in- 



tervals £ 



20 



< 



10 



> '^f^eff, the expressions ((^ . 



((2111 provide a good fit to the exact solution of equation 
((Tljl . The dirty limit formulas ((T3)l . ((TC|l are valid up to 
^20^ < 2^7// ^10 — '^^eff Figure 3 gives the ration of 
S,2/£i as a function of ^^"^ for a set oi ieff / £10 ■ At H ^ 
the oscillation length ^2 goes to infinity. Therefore the 
ratio is diverges at 0. With H increase the ra- 

tion rapidly decreases approaching the law C2/C1 °^ '^20^ 
at ^20 — 2- 




eff ^20 



FIG. 1: The decay length as a function of £20 calcu- 
lated for different values of . The open circles are the 
asymptotic curves, which have been calculated from 1201 for 
=2, 2.5 and 3. The open triangles are the asymptotic 
curves calculated from HSU for =1-1, 1-3, 1.5 and 2. The 
thin solid lines are the the asymptotic dependencies following 
from Eq. (1151 without the correction in the square brackets. 
These curves are calculated for =1.1, 1.3 and 1.5. 



C20 



coth 5 cos ' ^ 



^eff 



10 



^fo'sinhf^ 



_^ , arctan^ j (23) 
eff C20 / 



b = 



^20+^1 



^-isinhpr 



The discovered behavior of ^2 and ^1 is quite general 
and must be also observed in structures without ferro- 
magnetic ordering. An example is a normal filament of 
finite length, which is placed between superconducting 
banks and is biased by a dc supercurrent. It was shown 
[28j |. that the minigap induced to this filament from the 
S electrodes is not a monotonous function of phase dif- 
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FIG. 2: The difference between decay length and ^ 
a function of ^20 calculated for different values of shown 
in figure. The open circles are the asymptotic curves, which 
have been calculated from ( I2H for —2, 2.5 and 3. The 
open triangles are the asymptotic curves calculated from (I16II 
for =1.1, 1.3, 1.5 and 2. The thin solid lines are the the 
asymptotic dependencies following from Eq. llfcill without the 
correction in the square brackets. These curves are calculated 
for for =1.1, 1.3 and 1.5. 




FIG. 3: The ratio of oscillation and decay lengths ^2/^1 as a 
function of ^20 calculated for different values of ^^j^^. 



ference across the structure. This behavior could be also 
explained in terms of specific dependencies of ^2 and ^1 
upon electron mean free path in current biased systems. 



In summary, by solving the linearized Eilenberger 
equations we have calculated the real, ^1, and the imagi- 
nary, ^2 , part of a decay length as a function of exchange 
energy H and the mean free paths £, i^o, and for or- 
dinary, spin-orbit and spin-flip electronic scattering in a 
ferromagnet. These parameters, ^1, and ^2, characterize 
penetration of superconducting correlations into a ferro- 
magnet due to proximity effect and determine the decay 
and oscillation lengths of a critical current in long SFS 
Josephson structures. We have found the range of valid- 
ity of the expressions CO , @ , which are commonly used 
for interpretation of experimental data. In particular, 
the dirty limit expressions Q are valid if ^20^ < 0.5£~^j: 

fo'' ^10 — ^■'^^eff- The corrected expressions (fT5|) . (fTB|) 
can be used in a broader range of £,20 — '^^eff 
$10^ < '^^eff Further increase of exchange field makes 
the the length ^2 smaller than £e// thus breaking down 
the validity of approximations used in derivation of the 
Usadel equations. It is interesting to note that in certain 
parameter range the jumps occur in the dependence of 
S.2 vs ^20; while ^1 remains a continuous function of ^20- 



We have also demonstrated that the intuitive knowl- 
edge about the relation between ^1 and ^2, based on the 
dirty limit theory, has very limited range of applicabil- 
ity and can not be used for S^h > 5^ or for Ht > 0.1. 
In particular, an increase of H is not always accompa- 
nied by a decrease of ^1 and in a certain parameter range 
^1 may even increase with H. The fact that one may 
combine reasonably large decay length with the smaller 
period of oscillations looks rather attractive for possible 
applications of SFS Josephson junctions. 
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